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Abstract. To make predictions for the existence of "dark galaxies" , it is necessary to understand 
what determines whether a gas cloud will form stars. Star formation thresholds are generally 
explained in terms of the Toomre criterion for gravitational instability. I contrast this theory 
with the thermo-gravitational instability hypothesis of Schaye (2004), in which star formation 
is triggered by the formation of a cold gas phase and which predicts a nearly constant surface 
density threshold. I argue that although the Toomre analysis is useful for the global stability 
of disc galaxies, it relies on assumptions that break down in the outer regions, where star 
formation thresholds are observed. The thermo-gravitational instability hypothesis can account 
for a number of observed phenomena, some of which were thought to be unrelated to star 
formation thresholds. 
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1. Introduction 

Do all gas clouds with M ^> Mq contain stars? And if not, why do some clouds form 
stars whereas other do not? Observationally, it is nearly impossible to prove that a gas 
cloud has no stars at all, because a small number of old stars would generally be far 
too faint to detect. While it is difficult to put upper limits on the amount stellar mass, 
it is relatively easy to obtain interesting upper limits on the amount of star formation. 
This can for example be done through the non-detection of recombination lines such as 
Ha, which are generated in the H II regions around young stars, or through infrared 
observations. 

The H I components of most spiral galaxies are observed to extend beyond the region 
of active star formation. It has been known for some time, that star formation is generally 
limited to regions where the total gas surface density exceeds S c ~ 10 M Q pc~ 2 (e.g., 
Guidcrdoni 1987; Skillman 1987, Taylor et al. 1994). The existence of a surface density 
threshold for star formation is, however, usually explained in terms of the Toomre cri- 
terion for gravitational instability (e.g., Quirk 1972; Kennicutt 1989). Neither rotation 
nor pressure can stabilize a thin, gaseous, differentially rotating disk if the Toomre Q 
parameter, 

is less than unity (Safronov 1960; Toomre 1964); where the effective sound speed c s , the 
epicycle frequency ft, and the gas surface density S g all depend on radius r. According 
to this theory, star formation is suppressed when the gas surface density falls below the 
critical value 

£ c EE^Q£ g (1.2) 

7T(_r 

which depends on the local values of the velocity dispersion and the epicycle frequency. 

Kennicutt (1989) and Martin & Kennicutt (2001, hereafter MK01) tested the hypoth- 
esis that the Toomre criterion is responsible for the observed cutoff in the star formation 
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rate by measuring the ratio a = l/Q of the gas surface density S g to the critical surface 
density. Assuming a constant velocity dispersion of 6 kms -1 and after azimuthal aver- 
aging, MK01 found that a = S g /S c « 0.^3 at the truncation radius in a sample of 32 
nearby spiral galaxies. 

Elmegreen & Parravano (1994) emphasized that the need for a cold gas phase imposes 
a pressure threshold. Schaye (2004, hereafter S04) showed that the transition from the 
warm to the cold gas phase triggers gravitational instability on a wide range of scales. 
He demonstrated that this thermo-gravitational instability can account for the observed 
star formation threshold, which he predicted to fall in the range E c ~ 3-10 Mq pc -2 
(corresponding to a critical pressure P c / k ~ 10 2 -10 3 cm -3 K and a critical volume density 
7iH ~ 10 _2 -10 _1 cm -3 ), with weak dependencies on the metallicity, UV radiation field, 
turbulent pressure, and the mass fraction in gas. 

After briefly summarizing the S04 model in section [2l I argue in section [3] that the 
Toomrc Q hypothesis relies on assumptions that break down in the outer parts of galax- 
ies, where star formation thresholds are observed. In section [4] I discuss a number of 
observational findings that are naturally explained by the thermo-gravitational instabil- 
ity hypothesis. 

In the second part of my talk, I discussed the relation between star formation laws 
expressed in terms of surface densities, volume densities and pressures. I presented a 
method based on this theory with which it is possible to reproduce arbitrary Kennicutt- 
Schmidt laws using numerical simulations, without any free parameters. I do not have the 
space here to describe this model and its implications for the interpretation of simulations. 
For more on this topic, stay tuned for the paper "On the relation between the Schmidt 
and Kennicutt-Schmidt star formation laws and its implication for numerical simulations" 
(Schaye & Dalla Vecchia, in preparation). 



2. Thermo-gravitational instability 

The insterstellar medium of galaxies is highly complex. It contains a number of gas 
phases, including a cold (T <, 10 2 K), warm (T - 10 4 K), and hot (T > 10 4 K) 
phase. Feedback from star formation drives strong, supersonic turbulence which can both 
promote and suppress gravitational stability. Because of this complexity, it is difficult 
to predict exactly where and when star formation occurs. However, the outer parts of 
galaxies are much simpler. Here the gas is almost exclusively warm and there is no, or 
very little, star formation and hence turbulence is weak and typically subsonic. Moving 
in towards the center, we may therefore be able to predict where star formation begins, 
even if we cannot predict what happens further in. 

S04 noted that the UV background radiation implies the existence of a threshold surface 
density for the formation of a cold phase, which agrees with the empirically determined 
star formation threshold. He hypothesized that the drop in the pressure associated with 
the transition to a cold phase triggers gravitational instability and thus star formation. 
Conversely, for densities below the threshold, self-gravitating gas clouds are kept warm 
and stable by the UV radiation. 

To test this hypothesis, S04 constructed two model galaxies consisting of exponential 
discs embedded in cold dark matter halos. The discs were illuminated by UV radiation 
and contained metals and dust. For a fixed surface density, the temperature and volume 
density were determined numerically under the condition of vertical hydrostatic equilib- 

f We corrected their value of a, see Schaye (2004). 




Figure 1. Circular velocity as a function of radius for models HSB (left) and LSB (right). 
The short (long) dashed curve shows the contributions of the halo (disc). Also shown are the 
epicycle frequency (dot-dashed curves) and tt times Oort's constant A (dotted curves), both in 
kms -1 kpc -1 . Figure taken from S04. 




Figure 2. The Toomre Q parameter (solid curves), the temperature logT (dashed curves), 
the molecular fraction 5 + log /h 2 (dot-dashed curves) and the surface density log E(Mq pc~ ) 
(dotted curves) are all plotted as a function of radius for models HSB (left) and LSB (right). 
Figure taken from S04. 



rium. The radiative transfer and chemistry were calculated using the publicly available 
package CLOUDY (Ferland 2000). 

The parameters of the model galaxies were chosen to mimic a high and a low surface 
brightness galaxy. Figure Q] shows the rotation curves (solid), including the contributions 
of the halo (short dashed) and disc (long dashed), the epicycle frequency (dot-dashed) 
and tt times Oort's constant (dotted). 

The solid curves in Figure [5] show the Toomre Q parameter as a function of radius. 
The sudden drop in the Q values coincides with (and is caused by) a similar drop in the 
temperature (dashed curves), and a sharp increase in the molecular fraction (dot-dashed 
curves) . The transition to the cold phase, which coincides with the onset of gravitational 
instability (Q < 1), occurs at a fixed surface density. Note that the models become 
unrealistic shortwards of the critical radius, where feedback from star formation will 
increase the UV field, the mctallicity, and the turbulent pressure. 

Figure [3] again shows the Q parameter (thick solid curves) as a function of radius, but 
note that this time the y-axis is logarithmic. The value Q = 1 is indicated with solid, 
horizontal lines. A Toomre Q of unity implies marginal instability only for fluctuations 
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Figure 3. The Toomre Q parameter as a function of radius (thick solid curves) for models HSB 
(left) and LSB (right). The solid horizontal lines indicate Q — 1. The dot-dashed curves indicate 
the Q values below which the disc is unstable to perturbations with wavelengths of (top) 10 3 , 
(middle) 10 2 , and (bottom) 10 pc. The thick dashed curves show the Toomre Q parameter 
computed under the assumption that a — 6 kms -1 . The horizontal dashed lines indicate the 
critical Q value measured by Kennicutt (1989) under the assumption that a — 6 kms - . Figure 
taken from S04. 



with A = A C rit/2, where 



A« , (2.1) 



In general, instability of perturbations on scale A requires Q < Q C (A), where 



Q c (X)=2d- — . (2.2) 

V <^crit \ -Acrit / 

The dot-dashed curves in Figure [3] indicate the Q values below which the disc is unstable 
to perturbations with wavelengths of (top) 10 3 , (middle) 10 2 , and (bottom) 10 pc. The 
transition to the cold phase causes Q to drop far below unity, thereby triggering instability 
on a wide range of scales. 

Assuming a constant velocity dispersion of a = 6 kms -1 (thick dashed curves) would 
result in large errors in the Q-parameter. However, if we recalibrate the critical Q param- 
eter and compare with the value measured by Kennicutt (1989) (horizontal dashed lines), 
then we obtain agreement between the predicted critical radii (the intersections of the 
thick solid curves and the horizontal solid lines roughly coincide with the intersections 
of the thick dashed curves and the horizontal dashed lines). 

S04 computed a large range of models and found that the following function can 
describe the variation of the critical surface density with the parameters of the model 

log (W^zsl « 0.8+0.3 log(/) -0.3 log f-^-]+0.2 log (-—J——), (2.3) 



lM pc -2 / ' ' bKJJ ' b \Q.lZ Q J ' b \10 6 cm- 2 s^ 

where Z is the metallicity, / is the flux of ionizing radiation, and / = / s //thi with / g 
the mass fraction in gas and /th the fraction of the pressure that is thermal. The above 
formula (see S04 for second order corrections) was found to reproduce the results for a 
wide variety of galaxy models. 
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3. The assumptions underlying the Toomre hypothesis are invalid in 
outer discs 

In this section I will argue that the Toomre analysis is generally not applicable to 
star formation thresholds in the outer parts of disc galaxies, because the underlying 
assumptions are invalid. Below I will discuss the offending assumptions in turn. 

A Toomre Q smaller than unity is sufficient for star formation A Toomre Q 
of unity implies marginal instability only for fluctuations with A = A c rit/2. In the outer 
parts of disc galaxies, this critical wavelength is typically > 1 kpc (see Fig. [3]). Hence, 
while Q = 1 is relevant for large-scale modes such as spiral arms, we require fluctuations 
with much smaller wavelengths to become unstable to account for the formation of star 
clusters. Note that if E g ~ £ c (c = 6 kms" 1 ) is observed, then the condition Q C 1 
implies E c <C S c (cr = 6 kms -1 ) and thus a <C 6 kms -1 . 

The disc is thin In the outer parts of galaxies, the disc scale height is similar to the 
most unstable wavelength, thereby invalidating the assumption that the disc is thin. 

The perturbations are axisymmetric Beyond the truncation radius identified in 
azimuthally averaged observations, star formation does occur in isolated regions whenever 
the local gas surface density exceeds the global threshold (e.g., Ferguson et al. 1998). 
The Toomre analysis is valid only for axisymmetric perturbations and can therefore not 
straightforwardly be applied to such regions. 

The velocity dispersion is constant Although the Toomre analysis does not re- 
quire it, when it is applied to star formation thresholds it is generally assumed that 
the velocity dispersion is constant. Following Kcnnicutt (1989), most work assumes 
a = 6 kms -1 . For sufficiently large scale heights, the gas will be kept warm by the 
UV background radiation. Therefore, there will always be a component with a velocity 
dispersion of at least that of the warm phase. However, if there is any cold gas then the 
relevant velocity dispersion, and thus Q, can locally be much smaller. 

Rotation needs to be taken into account As mentioned above, to form stars, or 
even star clusters, we require instability for A <C A cr jt and hence Q <C 1. However, for A -C 
A C rit the Toomre criterion reduces to the two-dimensional Jeans criterion: A > Cg/G£ g . 
Hence, rotation is generally not directly relevant for star formation thresholds. 
The exception is the central part of the disc, where A cr ;t oc k~ 2 becomes small (see 
Fig. [1]) and the Coriolis force can, in principle, prevent collapse on scales relevant for 
star formation. However, in practice the disc cannot be considered thin in the central 
parts because the scale height is no longer small compared with the radius. Moreover, 
the bulge will need to be taken into account. 

4. Observations that are naturally explained by the 
thermo-gravitational instability 

A number of observational findings, some of which were thought to be unrelated, can 
be explained by the thermo-gravitational instability theory of star formation thresholds. 
Below I will discuss the most important ones. 

4.1. Why is the critical surface density constant? 

Why can a single constant, surface density threshold account for the observed star for- 
mation thresholds (e.g., Guiderdoni 1987; Skillman 1987, Taylor et al. 1994)? A constant 
threshold works works even on a pixel-by-pixel basis for high-resolution observations of 
nearby galaxies (de Blok & Walter 2006) and it can also account for the isolated star 
formation events beyond the optical radii of galaxies (e.g., Ferguson et al. 1998). 
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If the Toomre hypothesis holds, then S c oc c s n and should therefore vary from place to 
place. Even if we fix c s , as is often done, it should still depend on the epicycle frequency. 
If the thermo-gravitational instability is responsible, then the critical surface density will 
be nearly constant [see equation (|2.3p ] and the predicted value is in excellent agreement 
with the observations. 

Note that although S c is predicted to depend only very weakly on environmental 
factors such as the metallicity and the UV intensity, the differences may still become 
noticeable in extreme environments. For example, we predict S c 3> 10 M Q pc -2 for the 
first generation of stars at very high redshift (because Z <C O.IZq) and also very near 
strong sources of UV radiation. 

4.2. Why is a large fraction of disc galaxies inconsistent with the Toomre hypothesis? 

Why does a constant threshold work even when the Toomre criterion does not? Many 
studies have found so-called "subcritical discs". That is, galaxies for which Q(a = 
6 kms -1 ) exceeds the critical value throughout the region of active star formation (e.g., 
MK01; Wong & Blitz 2002; Auld et al. 2006). About 1 in 4 of the galaxies from the 
sample of MK01 fall in this category. As can be seen from the figures in these papers, 
the subcritical discs do agree with the constant star formation threshold hypothesis. 

4.3. Why is the velocity dispersion in the outer disc 8 kms -1 ? 

Why is the velocity dispersion of the outer HI disk constant and why does it have the 
value 8 kms -1 (e.g., Lo et al. 1993; Meurer et al. 1996)? Neither magnetically induced 
(Sellwood & Balbus 1999) nor gravitationally induced (e.g., Wada et al. 2002) turbulence 
can account for such a high (and constant) value. 

According to the thermo-gravitational instability picture, the gas in the outer disk, 
at surface densities below the threshold, is kept warm by the UV background radiation. 
Although this should be uncontroversial, theoretical studies have rarely taken radiation 
into account. 

4.4. Why does the same S c also work for tidal arms? 

Maybhate et al. (2007) found that the standard, constant surface density threshold can 
also account for observations of star clusters in tidal arms. Tidal arms can hardly be 
considered part of thin, differentially rotating discs, and the Toomre hypothesis is there- 
fore even less applicable to these observations than to those in the outer parts of disc 
galaxies. However, the thermo-gravitational instability does not assume a thin, rotating 
disc and should thus hold equally well in tidal arms. 

4.5. Why is gas with S g < S c nearly always warm? 

Why does the contribution of the cold interstellar phase to the gas surface density become 
negligible beyond the edge of the star forming disk (Braun 1997)? This observation 
is in excellent agreement with the predictions of the thermo-gravitational instability 
hypothesis: the formation of a cold phase triggers star formation and the two therefore 
go hand in hand. For £ < £ c the cold gas fraction, and thus the star formation rate, 
plummets. 

In the constant velocity dispersion version of the Toomre hypothesis there is no 1-1 
relationship between the ability to form a cold phase and the ability to form stars. While 
instability of the warm phase will presumably also eventually trigger the formation of 
a cold phase, the converse does not need to be true. Therefore, there is no reason why 
there could not be lots of cold gas in regions where there is no star formation. 
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4.6. Why does the H I surface density saturate at Shi ~ S c 

Many observations indicate that while the HI surface density increases with decreas- 
ing galactocentric radius, it typically saturates around 10 M Q pc~ 2 (e.g., Wong & Blitz 
2002). At smaller radii the total gas surface density can be much higher, but for £ g 3> 
10 M Q pc -2 most of the gas is molecular. Is it a coincidence that the H I surface density 
saturates at the star formation threshold? If the Toomre hypothesis is correct, then it 
would be. Although star formation may destroy the H I, there will be galaxies for which 
the disc is Toomre stable for £ g » 10 M Q pc~ 2 . 

In Schaye (2001) I pointed out that there is a physical upper limit to the HI surface 
density of gas clouds because the gas will quickly go molecular beyond a critical density, 
which depends weakly on the metallicity and the UV radiation field. The transition to 
a cold, molecular phase both limits the H I surface density and triggers star formation. 
Thus, the agreement between the critical surface density for star formation and the HI 
saturation limit is a prediction of the thermo-gravitational instability picture. 

It has been argued that MKOl's finding that the molecular fraction varies wildly at the 
critical radius for star formation is evidence against the thermo-gravitational hypothesis. 
However, a large scatter in the molecular fraction is exactly what is expected here. 
The critical radius R c is defined (by eye) as the radius within which star formation is 
ubiquitous and beyond which it is sporadic. The prediction is that for R » R c the gas 
is mostly atomic, while for R <C R c it is mostly molecular (provided that E g (r) is well 
behaved so that £ g < S c for R > R c and S g > E c for R < R c ). Since i? ~ i? c is the 
transition regime, we expect both atomic and molecular regions. Azimuthal averaging 
will then result in a wide range of molecular fractions. 

4.7. Why is the empirically determined critical Q different for irregulars? 

MK01 find that star formation is suppressed when S g /S c > a — 0.5. On the other hand, 
Hunter, Elmcgreen, & Baker (1998) measured a smaller value of a for a sample of irreg- 
ulars. However, Hunter et al. (1998) used a velocity dispersion of 9 kms -1 because the 
standard value of 6 kms -1 was inconsistent with the observations. Had they assumed 
used 6 kms -1 , they would have inferred a value of a that agrees within la with the Ken- 
nicutt (1989) value. This again reflects the fact that a constant surface density threshold 
describes the observations very well. If a is increased, then a constant surface density 
threshold implies a corresponding reduction of a. 

5. Conclusions 

S04 argued that star formation thresholds are set by the necessity to form a cold gas 
phase. Below a critical surface density the gas is kept warm by UV radiation. Conversely, 
when the surface density exceeds the threshold, the drop in the gas pressure associated 
with the phase transition triggers gravitational instability. As we discussed above, this 
hypothesis is in excellent, quantitative agreement with observations of the outer parts 
of galaxies. It can account for the observed star formation thresholds and the observed 
velocity dispersion, it can explain the lack of cold gas beyond regions of active star 
formation, and it also explains the coincidence between the critical surface density and 
the HI saturation limit. 

Does the Toomre Q play any role at all? It almost certainly does play an important 
role. Although it is not applicable to the outer parts, the molecular disc is thin compared 
with A cr it and can thus be characterized using the Toomre parameter. If Q < 1 then 
the disc is unstable to large-scale axisymmetric perturbations, which will for example 
lead to the formation of spiral arms. The resulting increase in the surface density can 
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trigger local instabilities, and not necessarily only through the formation of a cold phase. 
However, on scales relevant to the formation of star clusters, the 3-D Jeans criterion 
describes the local stability better than the local Toomre Q, and much better than the 
constant velocity dispersion version of the Toomre criterion. 
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Discussion 

TAYLOR: Recent GALEX observations do not show evidence for star formation thresh- 
olds. How do you interpret those observations? 

Schaye: The GALEX observations are very interesting. At first sight they appear to 
contradict the Ha observations, which do show evidence for star formation thresholds. 
GALEX detects star light, whereas the Ha comes from HII regions. GALEX is sensitive 
to star formation up to 100 Myr after the stars formed, whereas the Ha becomes difficult 
to detect after 10 Myr. The latter observations therefore offer more direct tests of theories 
of star formation thresholds. It is difficult to use GALEX observations to relate star 
formation to local gas densities because after 100 Myr feedback from star formation may 
well have destroyed/removed all the local H I and H2. GALEX may probe the cumulative 
effect of the isolated star formation events that are seen in Ha in the outer discs. 

Taylor: Stellar discs may not be truncated, although many do show evidence for a 
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change in the exponential scale length in the outer parts. Can star formation thresholds 
explain this? 

Schaye: Star formation does occur in the far outer parts of galaxies, when the local gas 
density exceeds the threshold. I could imagine some kind of change in the stellar surface 
density profile around the point where the azimuthally averaged gas surface density 
equals the threshold. Within that critical radius the star formation rate will be mostly 
determined by the efficiency of star formation rather than by the threshold. On the other 
hand, beyond the critical radius it will be determined more by the frequency of density 
fluctuations strong enough to exceed the threshold. However, the stellar disc is built up 
over long time scales, during which the critical radius may vary. Moreover, stars may 
move radially after they formed and stars in the outer discs may also have been accreted. 
These effects complicate the relationship between star formation thresholds and stellar 
surface density profiles. 



